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Abstract 

Introduction. Devices for collecting and storing energy from the external environment are low-power sources of electric 
energy that are actively used. The autonomous devices for monitoring the damaged condition of various structures include 
them as well. The working element of these devices is a piezoelectric generator (PEG) — a converter of mechanical 
energy into electrical energy. The design of PEG is associated with the preliminary construction of their mathematical 
and computer models, with the help of which the calculation and optimization of structures is carried out. One of the ways 
to model and calculate PEG is to develop approximate calculation methods based on applied theories. The applied theories 
for calculating bending vibrations of multilayer piezoactive plates are known and previously developed in the literature. 
However, in the scientific literature there is not enough information about bending and shear vibrations as a tool for 
improving the efficiency of engineering calculations of the described structures. The objective of this work was to develop 
an applied method for calculating bending and shear vibrations of piezoceramic plates, including porous ones. 

Materials and Methods. Piezoceramics PZT-4, including porous ones, were used as the piezoactive material of the plate. 
When using porous ceramics, the rigidity of the structure decreased to a greater extent than the piezoelectric modules, 
which made it possible to obtain a more effective PEG under mechanical action. The mathematical formulation was 
carried out within the framework of the linear theory of electroelasticity with plate polarization in thickness. The sides of 
the plate were electrodated, the right side was fixed, and a smooth contact in the vertical wall was set on the left side. 
Steady-state vibrations of the plate were caused by pressure on the front surfaces of the plate or the difference in electrical 
potentials at the electrodes. To calculate the characteristics of PEG, the authors proposed an applied theory based on 
hypotheses about the distribution of characteristics of the stress-strain state and the electric field. 

Results. Transverse vibrations of a piezoceramic plate in the low-frequency region (below the first bending-shear 
resonance) were studied. Due to the fact that the mathematical formulation was considered within the framework of the 
linear theory of elasticity, the problem was divided into the sum of two. The first one took into account the mechanical 
effect: a distributed load and a transverse force at the left end acted on the front surfaces of the plate, and the potentials 
at the electrodes were zero. In the second task, there were no mechanical loads, but the potential difference was set at the 
electrodes. Based on hypotheses about the distribution of deformations, mechanical stresses and electric potential, both 
problems were reduced to a system of ordinary differential equations and boundary conditions. Comparison with the 
results of calculations by the finite element method in the ACELAN package showed the adequacy of the proposed applied 
theory in the low-frequency region. 

Discussion and Conclusion. Since the formulation of the problem was considered in the linear theory of electroelasticity, 
and the low-frequency region was studied, the work succeeded in dividing the problem of bending-shear vibrations of a 
porous piezoceramic plate into two: bending — with mechanical action at zero potentials, and shear — when setting the 
potential difference and zero mechanical action. The corresponding hypotheses about bending and shear were used. Two 
systems of ordinary differential equations and boundary conditions, which were solved analytically without the use 
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of “heavy” finite element packages, were constructed. To compare the results and confirm the adequacy of the proposed 
method, the finite element modeling of such tasks was carried out in a specialized ACELAN package. The comparison 
showed that the error in determining displacements and electric potential when using this approach, in the case of setting 
mechanical loads and potential differences, did not exceed 6%. The method developed in the paper can be applied in the 
design of piezoelectric generators for energy storage in the low-frequency region. 
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AHHOTalna 

Beedenue. YctpolictBa cOopa u HakOMJIeHHA JHEP M3 BHEeIMIHeM CpebI MpeCTaBIAIOT COOOM MasIOMOMIHbIe HCTOU- 
HHUKH 3JICKTPHYeCKOM IHEPrHu, KOTOPbIe AKTHBHO HCIOMb3YIOTCA, B TOM YHCIIe B ABTOHOMHBIX IpHOopax MOHHTOpHHra 
MOBPeK CHHOTO COCTOAHHA pasIM4YHbIX KOHCTpyKUMH. PaOouuM 3IeEMeHTOM 3THX YCTPOHCTB ABIIACTCA IIbe309IEKTpH4e- 
ckui renepatop (IISI’) — mpeoOpa3oBatemb MexaHH4eckol 9HeprHuu B WIeKTpuyeckyro. Konctpyuposaune IIDT caa- 
3aHO C IIpeABapHTesIbHbIM MOCTpOeCHHeM HX MaTeMaTHYeCKHX HM KOMIBIOTEPHbIX MOeJeH, C MOMOLI[bIO KOTOPbIX IPOH3- 
BOJHTCA pacueT HU ONTHMU3AIIMA KOHCTpyKuMi. OHUM v3 crocoboB MoyemupoBaHua u pacueta IIOT saBnaetca pa3pa- 
OoTKa MIpHOMWKeHHBIX METOJOB pacueTa Ha OCHOBe IPUKIaTHbIx Teopuii. B mHTepaType U3BeCTHBI HM paHee pa3spa0oTaHbl 
IIpHKaHble TeEOpH pacueTa H3rHOHBIX KOIeOAaHHM MHOFOCOMHBIX Mbe3OaKTHBHBIX IWiacTHH. OqHako HHpopMalnH 06 
W3PHOHO-CABHTOBBIX KOeOaHHAX, KAK HHCTPyMeHTe MOBbILIeHHA I:pPCKTHBHOCTH HHKCHEPHBIX PaCUeTOB OMMCaHHBIX 
KOHCTpyKUMH, B HaydHOH MTeparype HegocTaTouHo. Llembio HacTosUel paOoTs! ABIIAIach paspaboTKa IpHkKaqHoro 
MeTOJja pacueTa H3rHOHBIX MH CABHTOBBIX KoJeOaHHi be3OKepaMHUeCKHX IIaCTHH, B TOM YMCIIe MOPHCTEIX. 
Mamepuanoi u memoooi. B xauectTBe 1be30aKTMBHOrO MaTepHalia MWaCTHHbI MCHOb3yeTCA Wbe3s0Kepamuka PZT-4, B 
TOM 4ucIe Mopuctas. [pu ucnoub30BaHHH MOpHCTOM KepaMUKH %KCCTKOCTb KOHCTpYKUMH YMeHbIUWaeTCA B OObIUeH CTe- 
TICHH, UCM MIbe3OMOJYIIN, YTO MO3BOIACT MOMYYNTh Ooee sspexTHBHBIM ITST mpu Mexanwyeckom BosyelicTBuu. Mate- 
MaTH4eCKad MOCTaHOBKa OCYIIeCTBIeCHa B PaMKaX JIMHeCHHOM TeEOPHU ZICKTPOYNPYFOcTH Ip NOMAPH3alMH TWIAaCTHHbI 
110 ToMMHHe. boKOBbIe CTOPOHBI MJIaCTHHI IIEKTPOAMPOBaHbl, IpaBad CTOPOHa 3akpelieHa, a Ha JIeBOn 3ayaH riayKn 
KOHTAKT B BEPTHKaJIbHOM CTeHKe. YCTaHOBMBIIHeCA KOJICOaAHHA IIaCTHHBI BbI3IBAIOTCA aBJICHHeM Ha JIMI[eBbIe IOBEpX- 
HOCTH IJIlaCTHHbI WIM pa3sHOCTbIO 3ICKTPHYeCKHX MOTCHUMAIOB Ha WIeKTponax. [1a pacueta xapaxtepucTuK IID B pa- 
OoTe MpeaaraeTca MpHKayqHad TeopuA, OCHOBaHHad Ha PHMOTesax O pacnipeseeHHU XapakTepHCTHK Halips»KeHHO-]e- 
(POpMUpOBaHHOTO COCTOAHHA HM 3IEKTPHYECKOTO MOA. 

Pe3yivmamot ucciedoeanua. PaccMoTpeHbl MonmepedHble KoIeOaHHA Ibe30KepaMH4eCKON MJIaCTHHbI B HA3KOUaCTOTHOM 
oOsacTH (HYKe TepBoro H3rMOHO-CABHTOBOrO pe3zoHaHca). B cusy Toro, YTO MaTeMaTHUeCKad MOCTaHOBKa paccMOTpeHa 
B paMKax JIMHeMHOM Teopuu yimpyrocTnH, 3ayjaya pa3zyeuacb Ha cyMMy Byx. B nepBolt yaHTEIBaIOCb MexaHH4eckoe 
BO3elCTBHe: Ha JIMI[eBble MOBEPXHOCTH IIaCTHHbI JelcTByeT pactipeyeseHHad Harpy3ka MU MomepedHas cia Ha JIEBOM 
KOHIe, a MOTCHUMasIbI Ha IICKTposax paBHbl HyO. Bo BTopo 3ayave MexaHH4eckHe Harpy3KH OTCYTCTBOBAJIH, HO 3a- 
JjaBalacb pa3HOCTb MOTeHUMaOB Ha 9eKTpopax. Ha ocHoBe rumoTe3 0 paciipeyzereHuu AepopMalHi, MexaHHyeckux 
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HallpsoKeHHH UM 9IEKTpH4eCKOrO MOTeHMasa OOe 3aa4H ObIM CBECHBI K CHCTeMe OOBIKHOBeHHBIX JHPepeHuMasIbHBIX 
ypaBHeHuli U TpaHW4HBIX yCOBMH. CpaBHeHue C pe3yJIbTaTaMM pacueTOB MeTOJOM KOHC4UHBIX 3JIEMCHTOB B TlakeTe 
ACELAN noka3asIv ayjeKBaTHOCTb MpesOKeCHHOM NpHKIaqHON TeOpHu B HU3KOUAaCTOTHOM OOacTH. 

O6cystcoenue u 3aknroyuenue. TlockolbkKy NOCTaHOBKa 3afa4 pacCMaTpHBasaCcb B JIMHEMHOM TeEOpHU BICKTPOyMpyrocTu U 
W3y4alacb HH3KOYAaCTOTHAaA OOACTb, B padoTe yAaslocb 3ayauy OO H3THOHBIX HM CABHTOBbIX KOJIeOaHHAX IWIaCTHHBI 13 T10- 
PUHCTOH Mbe30KepaMUKH pa3z{eIMTb Ha Be: H3rMOHy!O — C MeXaHHYeCKHM BO3eHCTBHeM IIpH HyJIeBbIX MOTeHUMasax Hi 
CJBHTOBy!lO — pH 3alaHuv pa3HOCTH NOTeHIMasIOB U HYJIEBOM MeXaHH4eCKOM BO3qecCTBHH. Ucnomb30BaHbl COOTBET- 
CTBYIOIUHe TuMOTesbI OO H3rHOe HM CABHTe, MOCTPOCHEI Be CHCTeMbI OOBIKHOBCHHBIX WHPPepeHUMaIbHbIX ypaBHeHHi 
TpaHW4HBIX YCHOBHH, KOTOpble PeWIalOTCA aHAMIMTH4eCKH 6€3 MCMOb3OBAHHA “TAKEJIBIX) KOHCYHO-3JICEMCHTHBIX MAaKeTOB. 
Ja cpaBHeHHA pe3yJIbTATOB HM MOATBEpYKTCHUA aeKBATHOCTH peIOKeCHHOTO MeTOIa NPOBeeCHO KOHCUHO-3JIEMCHTHOe 
MOJesIMpOBaHHe TakHXx 3a/ja4 B CIelHasIH3HpOBaHHoM Makete ACELAN. 3rTo cpaBHeHue noka3asi0, YTO OMHOKa B OMpesfe- 
JICHHM CMeNIeHHH UM SJIeKTpH4eCKOrO NOTeHHMasa Mp MCHOJ30BaHHU 3TOFO NOAXOA, B Cilyuae 3ajaHuA MeXAaHHYeCKHX 
Harpy30K H pa3HOCTH MOTeHUMaIOB, He lipeBbiiuaeT 6 %. PaspaOoTaHHbiii B CTaTbe METO MOX%KET ObITb IPHMeHEH IIpH IIpo- 
eKTHPOBaHHH Ibe303JIEKTPHYeCKUX TeHepaTOPOB HAaKOIWICHHA IHEPIuH B HU3KOUACTOTHOM OOACTH. 


Konouesbie cJIOBa: YCTpolcTBO cOopa 3Hepruny, TIbe309JIEKTPHAeCKHH reHepaTop, MOpHcTaaA KepaMHKa, w3rH6 
TW1aCTHHBI, CBHI HWIacTHHbI, WpuKilaqHaat TCOpuAa 
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OTHOIMICHHEe K CTAaTbe. 
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Introduction. Piezoelectric generators (PEG) are used to convert mechanical energy into electrical energy, followed 
by its accumulation. One of the application areas of PEG is the creation of low-power autonomous renewable sources of 
electric energy. The working element of the PEG is a piezoceramic element of a certain shape. The shape and type of 
deformation of this element determine the piezomodule, which characterizes the conversion of mechanical deformation 
energy into electrical energy. Thus, piezomodule d33 is associated with tension-and-compression along the axis of 
polarization, d3; — with the same deformation in the transverse direction to this axis, d:; — with shear. The use of porous 
ceramics makes it possible to create more efficient PEG. This is due to the fact that the elastic modules of porous ceramics 
decrease significantly stronger with increasing porosity than piezomodules. Thus, under the same mechanical load, the 
deformation amplitude of porous ceramics is greater; therefore, the output electric potential is also greater. 

PEG calculation can be performed by the finite element method implemented in ANSYS, ACELAN, COMSOL, and 
others packages. For piezoelectric elements, one or two sizes of which are significantly smaller than others (plates, rods), 
applied calculation theories can be constructed based on hypotheses about the distribution of mechanical and electric 
fields. Without the use of “heavy” finite element packages, applied theories make it possible to model various devices 
based on piezoactive materials. Piezoelectric, piezomagnetic and composite piezomagnetoelectric materials are 
considered as such materials. The construction of these theories is based on the acceptance of hypotheses about the 
distribution of mechanical, electric and magnetic fields. These hypotheses are related to the vibration mode of elastic and 
piezoactive elements of PEG. The most common designs are active and semipassive bimorphs based on multilayer plates, 
polarized in thickness with electrodes on the front faces, performing transverse bending vibrations. A number of papers 
are devoted to the study of devices with shear deformation of piezoelectric elements. An electric model with piezoelectric 
defining equations of mode djs and a single-degree-of-freedom model were combined to describe the energy collection 
characteristics of a piezoelectric cantilever in a shear mode in operation [1]. The proposed model is used to simulate the 
frequency dependence of the output peak voltage and power. The results show a good agreement with the experiment and 
the finite element calculation in ANSYS. In [2], a piezoelectric shear mode energy converter was developed to use the 
energy of a pressurized water flow. It converts the energy of the flow into electrical energy through piezoelectric 
conversion with vibration of the piezoelectric film. A finite element model has been developed to estimate the generated 
voltage of a piezoelectric film, which is in good agreement with the conducted field experiment. A one-dimensional fully 
coupled beam vibration model based on Timoshenko-type hypotheses, which provides a single common basis for energy 
analysis in shear and bending modes, is presented in [3]. In [4], the effect of the inhomogeneity of the properties of the 
plate under shear and torsional vibrations of its central part was studied. In experimental work [5], a multilayer cylindrical 
piezoelectric shear actuator (MCPSA) operating in shear mode ds, was presented for precision actuation under high 
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mechanical load. The actuator was made of piezoelectric ceramic rings Pb(Zr,Ti)O3 (PZT-51), which were concentrically 
assembled together in an electrically parallel connection with alternately positive and negative polarization in the axial 
direction. In [6], metamaterial of identical elementary cells was created, an artificial prototype of a device with 
characteristic patterned electrodes and piezoceramic subunits arranged in a row was designed and manufactured, which, 
as proven, ideally generated the synthetic shear deformation of the face. At the same excitation voltage, there was an 
increase in the displacement of the shear type by more than an order of magnitude, compared to the previous volumetric 
elements in mode djs. In the static formulation in [7], the field of electromechanical coupling in the shear-bending mode 
for an annular piezoelectric plate was theoretically established. In accordance with the classical theory of elastic plates of 
small bending and piezoelectric defining equations, an analytical solution of the bending deformation of the piezoactuator 
under the action of an electric field and a concentrated or evenly distributed mechanical load was achieved. The 
mechanism of generating bending deformation was explained by axisymmetric shear deformation, which additionally 
caused bending deformation of one piezoelectric plate in the form of a ring. This mechanism differs significantly from 
the mechanism of piezoelectric bimorphic or unimorphic drives, which were previously reported. The design of the 
annular piezoactuator has been optimized. In [8], a one-dimensional model was used to construct a sensor response 
function based on shear resonators (quartz cuts) of a bulk acoustic wave, which are promising for in-line measurements 
of fluid viscosity, e.g., in industrial processes. In [9], using the finite element method, a piezoelectric flying height control 
converter was investigated using a shear model deformation. In [10], the theory of a functionally graded plate 
with four-unknown shear deformations was used to express the displacement component. The distribution of the electric 
potential was a linear function of thickness. The plate was under mechanical load and electrical voltage. The basic 
equations and boundary conditions were derived using the virtual work principle. The analysis of stresses and 
deformations from the design parameters was performed. The electromechanical analysis of the stability loss of a 
piezoelectric nanoplate under shear using a modified theory of paired stresses with different boundary conditions was 
studied in [11]. To take into account the electrical effects, an external electrical voltage was applied to the piezoelectric 
nanoplate. A simplified theory of the first-order shear deformation was used. The basic differential equations were 
obtained using the Hamilton principle and nonlinear Von Karman deformations. Finally, the results showed that the effect 
of external electrical voltage on the critical shear load arising on the piezoelectric nanoplate was insignificant. In [12], 
using a combination of two classical approaches to modeling the nonlinear behavior of piezoelectric materials, a 
piezoelectric shear drive for an atomic-force microscope was investigated. Specifically, the novelty of the proposed 
method was in the fact that it combined two sources of nonlinearity of the field-dependent model from Miiller and Zhang 
with the frequency-dependent model from Damjanovic. The numerical results obtained using the finite element method 
(FEM) were compared to the experiment. 

Vibrations in which there is shear in addition to bending are less studied in the scientific literature, 1.e., piezomodule 
d\s, 1s the “working” one, whose value decreases with increasing porosity, but to a lesser extent than elastic modules. The 
latter circumstance makes it possible to build an efficient energy conversion device. Therefore, the development of an 
applied theory of PEG calculation using porous piezoceramics based on simplified models without “heavy” finite element 
packages seems to be a highly relevant task. The objective of this work was to construct an applied calculation method 
for steady-state transverse vibrations in the low-frequency region of a porous piezoceramic plate characterized by both 
shear and bending. 

Materials and Methods. The PEG under study was a piezoceramic plate (length /, thickness /), polarized in thickness, 
cantilevered on the right side, the left side was attached to an inertial mass that performed vertical vibrations and was 
fixed in the horizontal direction. The electrodes were located on the sides of the plate; therefore, with a potential difference 
on them and no mechanical load, the principal deformation in the low-frequency region was shear. Vibrations whose 
frequency was less than the frequency of the first resonance were considered. 

Mathematical formulation of the problem 

The mathematical formulation of the problem is described by a system of differential equations [13] and the 
corresponding boundary conditions. 

Pptitagpju-V-o=f; V-D=0 
o=cj--(€+Byé)—e; -E E=-Vo (1) 
D+caD=e;--(e+Gué)+95-E =(Vu+Vu")/2 


When considering porous ceramics of connectivity 3-0, equations (1) use effective physical constants determined 
through the ACELAN-COMPOS package [14]. These effective properties based on representative volumes (Fig. 1) were 
obtained in [15]. They are presented in Table 1. 
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Fig. 1. Representative volumes in ACELAN-COMPOS package [14] 


Table 1 
Effective properties of porous ceramics 

% of porosity 0 10 20 30 40 50 60 70 80 
P, kg/m? 7,500 6,750 6,000 5,250 4,500 3,750 3,000 2,250 1,500 
ee 10 Nin? 13.90 11.56 9.25 6.85 5.05 3.34 2.07 1.26 0.68 
e, 10”, Nim 7.78 6.15 4.66 3.14 2.10 1.16 0.62 0.28 0.13 
cn , 10%, Nin? 7.43 5.82 4.25 2.82 1.70 1.06 0.52 0.24 0.10 
ee 10" Nia 11.50 9.53 7.23 5.42 3.91 2.72 1.63 0.91 0.47 
chet 10! N/m? 2.56 2.23 1.83 1.44 1.10 0.74 0.44 0.23 0.10 
eo Clin 15.10 13.38 11.37 9.59 7.68 5.93 3.93 2.30 1.25 
es", C/m? 5.20 | -4.23 —3.14 | -2.07 —1.32 —0.75 0.43 —0.21 —0.10 
a Cm? 12.70 10.96 8.96 6.91 5.00 3.30 1.95 1.00 0.44 
Ki, / 86 730 663 582 509 439 349 263 191 122 

ic SEM 635 567 492 413 345 270 199 130 75 


According to Table 1, the dependences of effective elastic modules and piezomodules on the percentage of porosity 
[15] are plotted. They are shown in Figure 2. 
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Fig. 2. Dependences of values on percentage of porosity [15]: a — elastic models; b — piezomodules [15] 
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In accordance with the results shown in Figure 2, elastic modulus £; decreases significantly faster than 


piezomodule ds. 


Building an applied theory 
The developed calculation method consists of solving two problems: in the first, bending under the effect of a 


mechanical load with zero potential difference is considered; in the second, the shear caused by a potential difference in 
the absence of a mechanical load is modeled. In both cases, the absence of charges on the front faces of the plate is taken 
into account. Under the action of a mechanical load and an electrical potential difference, the results of the two tasks add 


up due to the linearity of the task. 
The action of the mechanical load, the potentials on the electrodes are zero. Hypotheses of the Kirchhoff-Love type 


are accepted regarding the equality of normal stresses to zero, and for displacements: 
U3 =UZ> (x), uw --(fuz (*)}2 (2) 
dx 
The distribution of electric potential over the thickness is assumed to be quadratic and symmetrical: 
2 
=D (x)e(22-1]0 +, (x)e(14 22) +, os 4} (3) 


To take into account the boundary conditions at the ends of the plate (x = 0, J), an expression for the transverse force 


is obtained: 


8033233 + 8e33 8 


C13€33 ; re 
1 (cugas +e) 803333 + 8e33 8 


OQ + €3) 


3 
12 Cc Cc +e% 7 h? 
33 ( 338 33 33h 


1 1 8¢33233 +8e33 4 d 
(casgss +03 )h x (4) 


2 2 
€33 (cisessh — €33€31h ) 


C13 C13 
1 (cssg +e% )h? es1(cises3h" —c33e31h" ) 
+ Cl 
12 C33 (casgss +e3)h? 
] eis (cixes3h? ~¢33e31h* )h a 
2 7 UZL2 (x) : 
24 C33 233 + €33 dx 


Taking into account the equality of the normal component of the electric induction vector on the front faces (z = +h) 
to zero, the equations for the unknown deflection UZ2(x) and distribution of the electric potential ®2(x) have the form: 


833233 +8e33 8 


C13€33 7 © 

1 (cssg35 +e%) as 8033233 + 8e33 8 1 
31 

2 C33 (casgss +03 )h? he 


y) 


2 
1 1 8033233 + 8e33 = 3 ag d : ®, (x) n 
2 2 (casgas +03 )h? h dx 


(5) 


2 2: 
€33 (cressh —€33e31h ) 


C13 2 2 C13 2 4 
(cssg33 +e )h e31 (c1ses3h —¢33e31h ) j 
C33 (casgas + €33 ) h 


e15 c13e33h" —c33e31h* h 4 
: | . : ~UZ2 (x) |-W* phUZ2(x)— p(x) =0, 
24 C33 233 + €33 dx 
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=0. 


The potential difference is set, the mechanical load is zero. Independence of lateral displacement from thickness, 
equality of the longitudinal displacement to zero, and the quadratic distribution of the electric potential over the thickness 


are assumed (3): 
U3 = UZ) (x), m1 =0. (7) 


Expression for the transverse force: 
d d 
Oi =-eu{ Luz, (:) |t-a[ er (3) (8) 


Taking into account the equality of the normal component of the electric induction vector on the front faces (z = +h) 
to zero, the equations for the unknown deflection UZ,(x) and distribution of the electric potential B2(x) have the form: 


2 2 
-ea| Suz o)r-es| So. 9} ~ p(x) -W*phUZ2(x) =0, (9) 
dx dx 


d° d° 
eis a2 (x) el 52 Dr (j=. (10) 


Research Results. The results of calculations based on the proposed applied theories were compared to the 
calculations of piezoelectric element vibrations (/ = 0.1 m, 4 = 0.01 m) at a frequency equal to 100 s" by the finite element 


method in ACELAN [16]. 
In the first problem, defined by equations (5) and (6), when specifying a uniformly distributed load p(x) = 1000 Pa-m 


and with boundary conditions: 
“uz (0) = 0, O1 |r=0= 0, D2 (0) =D, (/) = 0, UZ, (/) = 0, “uz: (/) = 0, (11) 


the following results were obtained, shown in Figures 3, 4. 
The calculations showed that the error in determining the vertical displacement was 5.8%, and for the horizontal 


displacement, it was 1.2%. 
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a) b) 
Fig. 3. Vertical displacement in ACELAN in the first problem: a — distribution; b — graph on the upper boundary 
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Fig. 4. Horizontal displacement in ACELAN in the first problem: a — distribution; b — graph on the upper boundary 


In the second problem, defined by equations (8) and (9), when setting zero load p(x) =0, potential difference 
Vo = 100 W and with boundary conditions: 
O} |x=0= 0, Bz (0) = 0, ®2 (7) =Vo, UZ2 (1) = 0, (12) 


the following results were obtained, shown in Figures 5, 6. 
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Fig. 5. Vertical displacement in ACELAN in the second problem: a — distribution; b — graph on the upper boundary 
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Fig. 6. Electric potential in ACELAN in the second problem: a — distribution; 
b — dependence on the longitudinal coordinate in the middle of the thickness 


The calculations showed that the error in determining the vertical displacement was 0.8%, and for the electric 
potential, it was less than 1%. It should be noted that the values of the horizontal displacement calculated in ACELAN 
turned out to be three orders of magnitude less than the maximum vertical displacement, which indicated the adequacy 
of the hypothesis (7). 

When setting mechanical loads and potential differences, the error of the proposed method turned out to be about 6 % 
for the displacement and electrical potential components. 

Discussion and Conclusion. As noted in the cited literature, the simultaneous use of bending and shear of a piezoelectric 
element can significantly increase its efficiency. In addition, the use of porous ceramics, due to different dependences of elastic 
modules and piezomodules on the percentage of porosity, also improves the output characteristics of PEG. 

In this paper, due to the linear formulation in the theory of electroelasticity, it was possible to develop an applied 
theory for calculating bending-shear vibrations of a piezoelectric element in the low-frequency region, which consisted 
of solving two problems: in the first, mechanical loads operated at zero potentials, and in the second, on the contrary, 
there were zero mechanical loads, and a potential difference was set. Based on various hypotheses about the distribution 
of mechanical and electric fields, two boundary value problems for systems of ordinary differential equations were 
obtained, which were solved analytically. The results of the calculation of displacements and electric potential were 
compared using the proposed method and the FEM implemented in the ACELAN package. These calculations validated 
the applicability of the proposed method, the error for which in calculating the above characteristics was 6%. This 
accuracy is sufficient for engineering calculations; therefore, the proposed method can be applied in the design of 
piezoelectric devices, including the collection and storage of energy. Further development of this applied theory will be 
aimed at covering a wider frequency range, including the first flexural shear resonance. 
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